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1.  INTRODUCTION 

I 

The  most  commonly  used  formulas  for  the  analysis  of  fiber-reinforced  and 
laminated  beams  are  those  given  by  the  familiar  "law  of  mixtures"  -£l^^In__ 
their  most  general  form,  these  formulas  are  derived  by  the  methods  of 
elementary  strength-of-materlals  theory,  on  the  basis  of  the  Bernoulli -Euler 
assumptions.  In  general,  of  course,  these  formulas  are  only  approximate:  for 
example,  when  calculating  the  deflection  of  the  simplest  type  of  laminated 
beam,  namely  a  sandwich.  It  Is  known  that  the  Inclusion  of  the  effect  of  shear 
Is  essential  cause  of  their  simplicity,  however.  It  Is  desirable  to 

use  them  whenever  one  may  do  so  with  sufficient  accuracy.  It  Is  therefore  the 
purpose  of  the  present  work  to  examine  the  validity  of  the  simple  formulas, 
and  to  determine  how,  and  when,  they  should  be  corrected.*.  More  specifically. 
It  will  be  shown  that  the  generalized  law  of  mixtures  Is  exact  for  the  case  of 
uniform,  uni -directionally  reinforced  beams  under  certain  simple  spanwlse 
bending  moment  distributions,  and  a  means  of  determining  the  required 
corrections  when  these  conditions  are  not  precisely  met  will  be  developed. 

The  work  presented  here  Is  the  counterpart,  for  reinforced  beams,  of  that  of 
[3].  Analogous  developments  have  been  carried  out  to  show  that  the  elementary 
formulas  were  valid  If  the  derivations  from  uniform  depth  are  small  and  smooth 
along  the  span,  respectively  In  [4]  for  homogeneous  and  In  [5]  for  fiber- 
reinforced  beams.  The  basic  method  In  these  works  was  first  Introduced  In  [6] 
to  study  the  effect  of  spanwlse  variations  of  temperature  distributions  on  the 
validity  of  the  elementary  formulas  for  the  case  of  rectangular,  beams,  and  was 
extended  In  [7,8]  to  beams  of  arbitrary  cross-section. 

In  the  analysis  presented  here  the  Inhomogeneous  beam  Is  considered,  at  the 
outset,  as  one  composed  of  a  single  material  but  with  variable  properties. 

The  basic  formulation  of  the  problem  thus  contains  spatial  partial  derivatives 
of  the  elastic  constants,  which  are  meaningless  at  the  Interfaces  between 
adjacent  layers,  since  sudden  jumps  In  the  material  properties  occur  there. 

As  one  Important  conclusion  of  the  present  work.  It  will  however  be  shown  that 
the  final  results,  both  for  stress  and  for  deformations,  can  be  expressed 
solely  In  terms  of  Integrals,  rather  than  derivatives,  of  the  moduli:  the 
difficulties  which  might  arise  because  of  discontinuities  In  material 
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properties  have  therefore  been  eliminated,  and  the  results  obtained  are  valid 
and  can  be  directly  used  for  reinforced  beams. 

It  will  be  seen  that  the  present  theory  leads  rather  quickly  to  fairly 
cumbersome  results,  particularly  If  the  geometry  of  the  beam  Is  at  all 
complex.  For  this  reason.  It  would  be  desirable  to  carry  out.  In  subsequent 
research,  a  certain  number  of  crucial  examples,  and  thus  to  be  able  to  reach 
some  hopefully  simple  criteria  for  the  Identification  of  cases  when 
corrections  to  the  elementary  formulas  are  Indeed  necessary.  Furthermore,  the 
development  of  a  computer  program  for  the  numerical  calculation  of  stresses 
and  deformations  from  the  techniques  developed  here  would  be  useful  as  well. 
Lacking  this  additional  work.  It  Is  felt  that  the  results  which  follow  present 
a  valuable  step  towards  the  accurate  and  practical  analysis  of  laminated  beams 
and  all  circumstances.  Some  further  details  on  the  present  work,  and  In 
particular  on  the  numerical  results,  may  be  found  In  [9]. 


2.  ELEMENTARY  RESULTS 


We  consider  In  this  paper  a  rectangular  beam  (occupying  the  space  0<x<L, 
-c<y<c,  0<z<w)  uni directionally  reinforced,  l.e.,  with  E  *  E(y)  and 
v>V(y).  The  elementary  formula  for  the  axial  stress  Is  the  generalized  law- 
of-mlxtures,  l.e.. 

The  corresponding  axial  displacement  u  Is  easily  calculated  from  the  equation 

u  -  4*  .while  the  curvature  Is 

•  ' 


where 


T  Is  the  temperature,  dA  ■  wdy  and  the  other  symbols  have  the  usual 
meanings.  Obvious  simplifications  occur,  for  example.  If  the  origin  of 
coordinates  Is  chosen  at  the  "centroid,"  l.e..  If  $  EydA  •  0. 


3.  BASIC  THEORY 


Under  the  normal  thinness  assumptions  (c/L«l»  w/L«l)  and  with  E  -  E(x,y) 
and  V-Y(x,y),  the  problem  requires  [6]  the  solution  for  the  Airy  stress 
functionary)  from  the  equation 

£  (1  2*  - +  .  1  &V  ?  21  (\**  \ 
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under  boundary  conditions  to  be  presently  discussed.  The  solution  may  be 
written  [4,5]  In  the  form 


If  the  quantities ^(x.y)  satisfy  the  equations: 


The  boundary  conditions  depend  on  the  details  ‘of  the  load  application.  We 
will  consider  for  simplicity  here  only  the  Isothermal  case  of  a  transverse 
load  q(x)  *  -d^M/dx^  applied  on  y  ■  c  (case  I  of  [3]);  the  boundary  conditions 
are  then: 


and  are  all  met  if  one  sets 


'**>(?,'■)  _  <L  .  ‘  _  c 


cvu  <j,  -  to 


(6) 


(6a) 

(6b) 


4.  STRESS  COMPONENTS 

Integration  of  eq.  (5a)  with  T  -  0  gives,  after  use  of  (6a)  and  of  Integration 
by  parts, 

^  v  i7) 

It  Is  evident  thatV^fl$y*  “^xxl  the  ^rst  term  In  a  series  similar 

to  (5)  for  the  axial  stress)  Is  Identical  with  the  elementary  one  of 
eq.  (1).  The  first  terms  of  similar  series  for  the  other  stress  components 
are  x 

V*  '^s;*  5*1*)  .(•*»]  m 

hMe»  4-  f'*J  m 

Further  terms  of  the  series  are  obtained  In  a  straightforward  manner  by 
Integration  of  (5b,c),  under  conditions  (6b),  Integration  by  parts  again  being 
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used  to  eliminate  all  partial  derivatives  of  the  moduli  with  respect  to  y. 

The  result  tends  to  be  rather  complicated:  for  example* Is  given  In 
Appendix  A,  and  Is  Indeed  expressed  solely  In  terms  of  moduli  and  their 
Integrals,  rather  than  derivatives  with  respect  to  y. 

For  the  Important  special  use  of  symmetrical  reinforcements  (l.e.,  with 
JyEdy  »<^yVdy  »  0,  such  as,  for  example,  an  ordinary  sandwich  beam)  the  stress 
components  take  the  form: 

A,*  -  -MrYr 

+C  -V  4  ... 

XV- 
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where  the  weighted  moment  of  Inertia  I*  -  w^Ey2dy  and  where 

-Mr^4 

v  j*;  h^^j) 

<v- v  hir^ 4 
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A  discussion  of  these  results  will  be  given  In  a  later  section. 

It  Is  Important  to  show  that  Is  Is  Indeed  possible  to  obtain  expressions  to 
all  values  of  1  which  are  free  of  derivatives  of  the  moduli  with  respect  to 
y.  To  prove  this  In  general  one  may  proceed  by  Induction,  l.e.,  by  showing 
that.  If  any  value  of  1  exists  for  which *f\„\  and  ft.g  and  their  first  two 
derivatives  with  respect  to  y  can  be  expressed  In  'a  form  which  Is  free  of  such 
derivatives,  then will  Itself  be  free  of  them.  This  follows  Immediately  by 
rewriting  (5c)  explicitly  forjj^&^and  noting  that  It  Is,  by  our  hypothesis, 
free  of  moduli  derivatives;  It  tl&n  follows  that*|&  andjpj  have  the  same 
character.  But  we  have  already  noted  that^  and<^2  likewise  have  this 
property,  and  thus  the  proof  Is  complete. 
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5.  DISPLACEMENTS 


The  strain  components  are  found,  once  the  stresses  are  known,  directly  from 
Hooke's  law;  from  these.  Integration  of  the  strain-displacement  relations  In 
the  usual  manner  [3]  gives  the  desired  displacement  components.  The  series 
for  the  axial  displacement  u  Is  of  the  form 


where  the  second  and  third  terms  (l.e.,  u2  and  u3)  are  given  expllctly  In 
Appendix  A.  The  expression  for  the  curvature ^v/^x*  is  as  follows: 


6.  DISCUSSION  OF  RESULTS 

Several  observations  regarding  the  results  obtained  above  may  be  made: 

1.  All  results  obtained  for  non -homogeneous  beams  contain  no 
spatial  derivations  of  the  moduli.  They  are  thus  In  the 
desired  form,  and  no  computational  difficulties  arise  at  the 
Interfaces  between  component  materials. 

2.  At  each  Interface  the  shear  stresses,  the  stresses  normal  to 
the  Interface  and  the  transverse  displacements  are  continuous 
(as  they  should  be),  but  the  axial  normal  stress  and  the  axial 
displacement  are  not. 

3.  The  results  are  all  In  series  form,  successive  terms  of  which 
are  proportional  to  progressively  higher  spanwlse  derivations 
of  the  bending  moment  and  on  higher  powers  of  the  thickness 
dimensions.  They  are  therefore  entirely  analogous  to  the 
corresponding  ones  for  homogeneous  beams,  and  hence  the 
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conclusions  well  known  for  that  case  can  be  Immediately  carried 
over  to  reinforced  beams.  In  particular,  note  that  the  law-of- 
mlxtures  (l.e.,  the  first  term)  Is  complete  solution  for  the 
axial  stress  and  for  the  curvature  If  the  bending  moment  Is  of 
the  form  M  ■  a  +  bx.  The  sum  of  the  first  and  second  terms  Is 
the  complete  solution  of  M  is  at  most  a  cubic  polynomial  In  x, 
and  In  general  the  accuracy  of  the  elementary  formulas 
Increases  as  one  considers  thinner  and  thinner  beams  and 
smoother  and  smoother  moment  distributions.  For  further 
discussions  of  the  meaning  of  "smooth"  and  of  further 
conclusions,  the  reader  Is  referred  to  the  previous  treatments 
(e.g.,  [3]  or  [7]).  Of  course,  all  the  present  results  reduce 
to  the  known  ones  valid  for  homogeneous  beams  If  simplified  by 
taking  both  E  and  V  to  be  constants. 

4.  In  the  calculation  of  stresses  In  sandwich  and  multi-layer 
beams,  the  terms  other  than  the  first  one,  as  can  be  expected, 
are  of  very  minor  Importance.  This  Is,  of  course,  not  true  In 
the  calculation  of  deflections,  where  the  second  term,  for 
example.  Includes  the  effect  normally  referred  to  as  shearing 
deformations  and  Is  therefore  of  substantial  magnitude  when  one 
of  the  component  materials  Is  weak  In  shear.  Such  matters  as 
the  bending  rigidity  of  sandwich  faces,  and  the  normal-stress 
carrying  capacity  of  the  sandwich  core  are  automatically 
Included  In  the  present  theory;  they  are  normally  (cf.  [2]) 
neglected  In  sandwich  analyses,  and  their  effect  can  therefore 
be  estimated  with  the  aid  of  the  present  formulas. 

5.  In  complete  analogy  again  with  homogeneous  beam  results,  the 
correction  terms  In  the  series  are  often  more  Important  for 
thermal -stress  analysis  than  for  Isothermal  beams  [7].  This 
occurs  because  non-trlvlal  cases  In  which  the  first  term 
vanishes  can  arise,  as  for  example,  that  of  a  free  multi-layer 
beam  under  a  temperature  linear  with  y.  If  the  various  layers 
have  different  moduli  but  the  same  coefficient  of  expansion. 

In  such  cases  of  course  the  entire  stress  Is  given  by  the 
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"correction"  terms,  and  they  cannot  therefore  be  neglected, 
although  no  detailed  consideration  of  the  thermal  case  Is 
Included  here. 


REFERENCES 


1.  Broutman,  L.  J.  and  Krock,  R.  H.,  Modern  Composite  Materials,  Addison* 
Wesley  Publishing  Co.,  Reading,  Massachusetts,  1^7. 

2.  Plantema.  F.  J..  Sandwich  Construction,  John  Wiley  &  Sons,  Inc.,  New  York, 
1966. 

3.  Boley,  B.  A.  and  Tollns,  I.  S.,  "On  the  Stresses  and  Deflection  of 
Rectangular  Beams,"  Journal  of  Applied  Mechanics,  23,  pp.  339,  1956. 

4.  Boley,  B.  A.,  "On  the  Accuracy  of  the  Bernoulli -Euler  Theory  for  Beams  of 
Variable  Section,"  Journal  of  Applied  Mechanics,  30,  p.  373,  1963. 

5.  Boley*  B.  A.  and  Testa,  R.  B«,  "Thermal  Stresses  In  Composite  Beams," 

International  Journal  of  Solids  and  Structures,  vol.  5,  1969,  pp.  1153- 
TO - 

6.  Boley,  B.  A.,  "The  Determination  of  Temperature,  Stresses,  and  Deflections 
In  Two-dimensional  Thermoelastic  Problems,"  Journal  of  Aeronautical 
Science,  23,  1956,  p.  67, 

7.  Boley,  B.  A.  and  Weiner,  J.  H.,  Theory  of  Thermal  Stresses,  John  Wiley  & 
Sons,  Inc.,  New  York,  1960. 

8.  Barrekette,  E.  S.,  "Thermoelastic  Stresses  In  Beams."  Journal  of  Applied 
Mechanics,  vol.  27,  Trans.  ASMS,  vol.  82,  Series  E,  1960,  pp.  465-4/3. 

9.  Doelng,  F.  M.,  "Stresses  and  Deflections  In  Layered  and  Sandwich  Beams," 
M.S.  Thesis,  Cornell  University,  Ithaca,  New  York,  May  1973. 


APPENDIX  A 


The  expression  for  «^2  for  the  case  1n  wh1ch  E  ■  E(y)  and  V«  v'(y)  is  obtained 
by  the  process  described  In  Section  3  to  be  as  follows: 


+  1 
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Here 


The  stress  components  obtained  from«p2  are  now  6as11y  calculated  (cf.  [9]), 
but  are  not  explicitly  listed  here  for  the  sake  of  brevity.  They  form  the 
first-order  correction  to  be  app7'ed  to  the  stress  components  listed  In 
eqs.  (7)  and  (8).  It  Is  easily  seen  that  all  these  stress  components  are 
themselves  free  of  derivatives  with  respect  to  y. 


Eq.  (Al)  and  the  corresponding  stress  components  reduce  Identically  to  those 
of  [3]  for  the  special  case  of  homogeneous  beams.  Simplifications  resulting 
In  the  case  of  symmetrically  placed  reinforcements  are  given  In  eq.  (9). 


The  displacement  components  are  obtained  by  means  of  the  steps  outlined  In  the 
text.  For  the  axial  components,  the  result  Is: 


